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ABSTRACT 



Some X-ray dim isolated neutron stars (XDINS) and central compact objects in supernova remnants (CCO) contain absorption features 
in their thermal soft X-ray spectra. It has been hypothesized that this absorption may relate to periodic peaks in free-free absorption 
opacities, caused by either Landau quantization of electron motion in magnetic fields B < 10" G or analogous quantization of ion 
motion in magnetic fields B > 10'^ G. Here, I review the physics behind cyclotron quantum harmonics in free-free photoabsorption, 
discuss different approximations for their calculation, and explain why the ion cyclotron harmonics (beyond the fundamental) cannot 
be observed. 
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1. Introduction 

Thermal radiation from neutron stars can provide important 
information about their physical properties. Among neutron 
stars with therma l-like radiation spectra (see, e.g., reviews by 
iKaspi et al.ll2006l and lZavlin,2009i) . there are two classes of ob- 
jects of_2^ticularinterest: central compact objects (CCOs; see, 
e.g., Ide Lucall2008l) in supernova remnants and X-ray dim iso- 
lated neut ron stars (XDI NSs, or the Magnificent Seven; see, e.g., 
review bv lTurollall2009l) . 

The CCOs are young, radio-quiet isolated neutron stars with 
relatively wea k magnetic fields B ~ (lO"^ - lO") G (e.g., 
iHalpern & Got thelf 2010 and references therein). The XDINSs 
are older and are believed to have much stronger fields B > 
(10'^ - 10'^) G (lHaberlll2007t Ivan Kerkwiik & Kaolanl |2007[ 
lTurollal2009l) . For some CCOs and XDINSs, there are estimates 
of B, and in some cases only upper limits to B are available. 

In the past decade, broad absorption lines have been detected 
in the thermal spectra of several i solated neutron stars (see, e.g., 
! Ivan Kerkwiikl2004tlHab erl 2007: Ivan Kerkwiik & KaDlanl2007[ 
' and references therein). In all but one case, the energies E'^ of 
the absorption are centered on the range 0.2-0.7 keV and the 
effective black-body temperatures are T^^^^ a; 0.1 keV. Here 
and hereafter, the Boltzmann constant is suppressed, and the su- 
perscript "oo" indicates a redshifted value. In particular, it has 
been found that: (i) the spectrum of RX J1605. 3-1-3249 with 
^eff bb ~ 96 eV has a broad absorption at E'^ ^ 0.4-0.5 keV and 
a possible second absorp tion at 0.55 keV (Ivan Kerkwiik et al] 
l2004l:lvan Kerkwiikll2004l) : (ii) RX J0720.4-3125 exhibits an ab- 
sorption feature at ^ 0.27 keV ( Haberl et al. 2004) and a 
possible second absorption at 0.57 keV (iHambarvan et al.l2009l) . 
while the effective black-body tem perature varies ov er years 
across the range T^^^^^ ^ 86-95 eV ( iHohle et al.|[2009l) : (iii) the 
spectrum of RBS1223 (RX J 1308. 6 + 2127) was reproduced by 



a model with T^^^^ ^ 102 ± 2 eV and two absorption lines at 
E^ ~ 0.3 keV and E^ ~ 0.6 keV dSchw ope et al. 2007); and 
(iv) the spectrum of RBS1774 (IRXS J2 14303.7 + 065419) with 
TZ: * 102 eV shows indications of a line at E^ ^ 0.3-0.4 keV 
and an absorption edge at 0.7 3 -0.75 keV jCropper et al1l2007l ; 
Kap lan & van KerkwiikI 120091: ISchwope et alJ l2009l) . The first 
discovered isolated neutron star with absorption lines, CCO IE 
1207.4-5209, has t wo absorption featur es centered on E'^ » 0.7 
keV and 1.4 keV (ISanwal etalJl2002h and an effective black- 
body temperature (which may be nonuniform) of T' 



eff.bb 



0.16- 

0.32 keV dZavUn et alJlT998l: Ide Luca et al. I I2OOI . For this ob- 
ject, two more harmonically spaced absorption fe atures (at E^ ^ 
2. 1 k e V and 2.8 keV) wer e tentatively detected (iBignami et al.l 
I2OO3I : Ide Luca et a l.''2004'). but were later shown to be statisti- 
cally insignificant (Mori et al. 2005). We note that realistic val- 
ues of the effective temperature T^, obtained using atm osphere 
models, can differ from by a factor < 2-3 (see, e.g.. lZavlinl 

2009, and references th erein). 

Many authors (e .g., [Sanwal et al.ll2002l:lBignami et aLllIool 
Ide Luca et al.ll2004l) have considered the theoretical possibility 
that the absorption lines in the thermal spectra of the CCOs and 
XDINSs may be produced by cyclotron harmonics, formed be- 
cause of quantum transitions between different Landau levels of 
charged particles in strong magnetic fields. IZane et alJ (1200 ih 
discussed this possibility prior to the observational discovery 
of these absorption features. The fundamental cyclotron energy 
equals 



ho)^ = heB/mc = 11.577 Bn keV 



(1) 



for the electrons and hn^ = hZeB/Mc = 6.35 (Z/A)B 12 eV for 
the ions, where m and M are the electron and ion masses, re- 
spectively, Z and A are the ion charge and mass numbers, and 
B12 = B/IO'^ G. In the following, we consider protons, whose 
cyclotron energy is 
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h£lc = heB/Mc = 6.305 Bn eV. 



(2) 
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Be ginning with the pioneering work of iGnedin & SvunvaevI 
(11974 . numerous papers have been devoted to the physics 
and modeling of cyclotron l ines in X-ray spec tra of 



1977 



1999; 



accreting neutron stars (e.g. , Daueher tv & Ventural 
Pavlov etal.l 1 19801 IWangetaLl 119931: L^ va & Hardin J 
Arava-Gochez & Harding| |200tf, 'Nishimura 2005 u2008i) . These 
emission lines have b een observed in man y works follow- 
ing their discovery by iTriimper^ al.i (Il978h . Cyclotron har- 
monics have bee n found in spectra of se v eral X-ray pulsars 
in binaries (e.g. lRodes-Rocaetal.1 12009L lEnoto et alJ l2008t 
iPottschmidt et al.l l2004. and references therein), and up to four 
harm onics were registered for one of them CSantangelo et alj 
Il999h . 

In the photospheres of isolated neutron stars, unlike X-ray bi- 
naries, the typical energies of charged particles are nonrelativis- 
tic. In this case, first-order cyclotron transitions of free charged 
particles are dipole-allowed only between neighboring equidis- 
tant Landau levels and form a single cyclotron resonance with no 
harmonics. Special relativity and non-dipole corrections at the 
energies of interest can be estimated to be max(reff, £'a)/'wc^ ~ 
10"^ for the electrons and max(reff, E.^)IMc^ < 10"^ for the pro- 
tons. 

Beyond the first order in interactions, transitions between 
distant Landau states are also allowed in the nonrelativistic the- 
ory. They are, in particular, caused by Coulomb interactions 
between plasma particles. Thus cyclotron harmonics appear in 
free-free (bremsstrahlung) cross-sections. To obtain £3 ~ 0.1- 
1 keV, one may assume either the electron cyclotron harmonics 
at B ~ 10'*'-10" G, according to Eq. ([T), or proton cyclotron 
harmonics aXB~ lO'^-lO'^ G, according to Eq. (|2]i. 

iPavlov & Shibanovl (1 19781) presented the calculations of 
spectra for isolated neutron stars with prominent electron cy- 
clotron harmonics due to t he free -free absorption in the atmo- 
sphere. ISuleiman ov et al.l (1201 Obi) performed a si milar atmo- 
sphere modeling and concluded, in agreement with IZane et all 
(|200l|), that electron cyclotron harmonics could be observed in 
ceo spectra. Proton cyclotron harmonics cannot be calculated 
based on the assumption of classical proton motion, used by 
these authors. 

In this paper, I review the physics of free-free photoabsorp- 
tion in strong magnetic fields, discuss restrictions on different 
published approximations for free-free opacities, and present nu- 
merical results that demonstrate the relative strengths of the elec- 
tron and proton cyclotron resonances under the conditions char- 
acteristic of the atmospheres of isolated neutron stars with strong 
magnetic fields. This gives a graphic explanation of the small- 
ness of the ion cyclotron harmonics. I also demonstrate that the 
contribution of bound-bound and bound-free transitions to the 
opacities of neutron stars with B > 10'^ G is much larger than 
that of the proton cyclotron harmonics. 

In Sect.|2] quantum mechanical integrals of motion and wave 
functions of a charged particle in a magnetic field are recalled 
for subsequent use. Section [3] is devoted to the properties of an 
electron-proton system in a magnetic field that is quantizing for 
both particles: general equations for calculation of wave func- 
tions are given, and the Born approximation is considered in de- 
tail. In the same order, general expressions and Born approxi- 
mation are considered in Sect.|4]for photoabsorption matrix el- 
ements and cross-sections. Section |5] gives numerical examples 
of cyclotron harmonics in free-free photoabsorption with discus- 
sion and comparison of various approximations. Consequences 
for the CCOs and XDINSs are discussed in Sect.|6] and Sect.|2] 
presents our summary. 



2. Charged particles in a magnetic field 

Since special relativity effects are of minor importance in the at- 
mospheres of isolated neutron stars, we use nonrelativistic quan- 
tum mechanics. 

We assume that the magnetic field vector B is along the z 
axis and consider its vector potential the cylindrical gauge to be 



A{r)^^-Bx{r- 



rA) 



(3) 



with an arbitrary center in the xy plane. 

We recall the description of a charged particle in a 
uniform magnetic field (e.g., 'Johnson & L ippmaiinI Il949t 
iLandau & Lifshitzlll976HJohnson et al...l983.) . The Hamiltonian 
equals the kinetic energy operator 



H 



(1) 



mr 

where m is the mass. 



mr = p- (Q/c) A(r) 



(4) 



(5) 



is the kinetic momentum, Q is the charge, and p is the canonical 
momentum conjugate to r. In Eq. (|4| and hereafter, "±" denotes 
the "transverse" part, related to the motion in the xy plane. 

A classical particle moves along a spiral around the normal 
to the xy plane at the guiding center re- In quantum mechanics, 
rc is an operator, related to the pseudomomentum operator 



hk - mr + (Q/c) Bxr, 



(6) 



where - (c/QB^)hk X B. Its cartesian coordinates (xc,yc) 
commute with //^'\ but do not commute with each other: 
[■^cjc] = -ihc/QB. Another important integral of motion is 
the z-projection of the angular momentum L = (QB/c) (r^/2 - 
H^l^ /mcdil), where a>c — \Q\B/mc is the cyclotron frequency. 

The eigenvalues of H^^^ are given by - (n + ^)hojc, where 
n = 0, 1, 2, ... is the Landau quantum number The simultaneous 
eigenvalues of f, are (sign Q) hs with integer s > -n, and eigen- 
values of the squared guiding center r^ equal to a^(2s + 2n + 1), 
where = (fic/\Q\By^^ is the so-called magnetic length. 

In general, H^l^ should be supplemented by (-B ■ jj), where 
fi - gmag{e/2mc)S is the intrinsic magnetic moment of the 
particle, 5 is the spin operator, and ^mag is the spin ^-factor 
(^mag = -2.0023 and 5.5857 for the electron and the proton, 
respectively). In most applications, one can choose the repre- 
sentation where the electron and proton spins have definite z- 
projections +h/2 and set the electron ^-factor to -2, thus re- 
garding the excited electron Landau levels as double degenerate. 

The form of a wave function depends on a choice of the 
gauge for A(r). We consider the cylindrical gauge given by 
Eq. (13) centered on the coordinate origin {rA = 0). The eigen- 
functions of and £r in the coordinate representation are 



%,.,AXr) = -T72 X 



<i>nAr±), if e < 
<D* (r J, if e > 



(7) 



where k, = p^jhis the particle wave number along the field, 
is the normalization length, r - {x, y) - (r^ cos 0, r^ sin 0), 



yl2^L 



(8) 



is Landau function, the asterisk denoting a complex conjugate , 
and In'niu) is a Laguerre function (e.g.. Sokolov & Ternovll98^ . 
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We define cyclic components of any vector a as a-ti - 
(Gx + ioy)/ V2 and aq = a^. The transverse cyclic components 
of the kinetic momentum operator given by Eq. (|5]l transform 
one Landau state |n, s)^, characterized by <l>J,*5'(rx), into another 
Landau state 



where 



= — —B(rxp) + HK 

2m, Mc ^ 



(16) 



tL I 

mra \n, s}_i_ — -yn + 1/2 + a/2 \n + a, s — a)j_, 

where a = +1, a = a if Q < 0, and a = -a if 2 > 0. 



(9) 



is the Hamiltonian of the harmonic motion in the xy plane, p is 
the momentum conjugate to r. 



n = p+ — Bxr, 

2c 

Kl 



(17) 



3. Electron-proton system in a magnetic field 

The Hamiltonian of the electron-proton pair (i.e., of H atom) is 



where rep = /"e - rp. The kinetic part can be written as 



(10) 



(11) 



Hk = — + — Kb-(Bx r), and Kg ^ K + - B x j(i8) 
2mji mnc c 

We note that Hk - when rg - r^. 

The first term in Eq. ( fTsT l is the total kinetic energy along 
z, uncoupled from the relative electron-proton motion, therefore 
we set P, = without loss of generality. 

The eigenvalues of equal E^.. = + /i^, where n > 
and = « + .?> are the electron and proton Landau num- 
bers, respectively, and hs are eigenvalues of the relative angular 
momentum projection operator ((.^p - /'.^). 

We construct numerical solutions of Eq. (fl4l i in the energy 
representation for rg - tjTc ij] e [0, 1]) in the form 



where m and M are the electron and ion masses, R - 

(M/mii)rp + (m/m}i)re is the center of mass, me - m + M is lA^V) = '^n's'ir±) g^^,l,,.^{z) 
the total mass, and m« - niM/mn is the reduced mass. n's' 

Since the electron and the proton have opposite charges, their 
orbiting in the transverse plane is accompanied by a drift across 
the magnetic field lines with velocity Vdrift depending on the dis- 
tance between their guiding centers or equivalently on the total 
pseudomomentum K 



(19) 



re,. 



BxK, K^P 



where k is the composite quantum number enumerating quan- 
tum states. One retains in Eq. (fT9T l as many terms (n = 
0,1,2,..., Hmax; - -n,—n+ 1,-n + 2, . . . , i^ax) as needed to 
reach the desired accuracy. We choose a principal ("leading") 
term (n, s) and define "longitudinal" energy of the state \k) as 

(';) 



(12) Ef. - Ei; - E^^. The functions gj,\,.i^(z) are computed from 



where P is the canonical momentum conjugate to R. 

In quantum mechanics, it is not only true that the pseu- 
domomentum operator K - hk^ + hkp commutes with H, 
but also that its cartesian components {Kx,Ky,K,) commute 
with each other. Therefore, all components of r^. can be 
determined simultaneously. Coordinate eigenfunctions of the 
pseudomome ntum operator with eigenvectors K are given by 
dGor'kov & b zyaloshinskii 1968) 



T- TT + yn"s",n"s"{rB,z) + W s" I Hk I n"s"}^ 



+Eb, 



E^ 



gn"s";K{z) 



^{r^,rp) = exp 



)^[K.'-Bxnp).R 



(13) 



From the general Schrodinger equation W¥ - E^, one can de- 
rive an equation for tfrKirep), which has the form 

H,,i,pK^Eif/K, (14) 
where the effective Hamiltonian H^ei depends on K. 

3.1. Exact solution 

Solutions of Eq. ( fT4l i for arbitrary K in strong magnetic fields 
in the cylindrica l gauge represented by Eq. (O were obtained 
by IVincke et all ([l992) and Potekhin (1994) for bound states, 
and by iPotekhin & Pavlov! (Il997h for continuum states of the 
electron-proton svstem. 'PotekhinI (1 19941) used the variable r - 
rep - rs as an independent argument of the wave function and 
found that the most convenient parametrization in Eq. Q is - 
[(M - m)/mu] rg. Then 



H 



rel 



-I- -^+H^ 



2wh 2m t, 



(15) 



-Y^[v„",',n'ArB,z) + {n"s" \Hk\ n' s')^)g„,,.Jz), (20) 

where n" - 0, 1, 2, . . . , n„iax, s" - -n",-n" + 1,-n" + 
2, ...,imax, 2^' denotes the sum over all pairs {n',s') except 
(«", s"), and 

Vn"s".n'ArB,z) = («"/' | " l\r + r^l | n's')^ (21) 

are effective potentials (see |Potekhinlll994l for calculation of 
these potentials and matrix elements {n"s" \ Hk \ n's'}±). 



3.1.1. Bound states 

Bound states of the H atom can be numbered as \k} = \nsvK}, 
where v = 0, 1,2, . . . enumerates energy levels for every fixed 
pair (n, s) and controls the z-parity according to the relation 
gn',s';f<i-z) = (-iygn',s';Az). The longitudinal energies eI = 
eI„(K) are determined from the system of equations in Eq. ( l20b 
together with the longitudinal wave functions gn's'Az). 

The atomic states are qualitatively different for small and 
large K± values. For small K^, the electron remains mostly 
around the proton, the energy dependence on K± is nearly 
quadratic, so that the transverse velocity Vdrift - dE/^/dK^ 
is nearly proportional to K±, i.e., Vdrift ~ K±/m^. The ef- 
fective mass exceeds ma and increases with increasing B 



4 
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(|Vmcke&Bave"l988"). For large K^, the atomic state is decen- 
tered (Gor'kov & Dzvaloshinskii 1968): the electron finds itself 
mostly around r^, rather than around the proton. In the latter 
case, Vdiift decreases with increasing K^. The two families of 
states are separated by the critical value of the pseudomomen- 
tum, /Tcr ~ (2mH|-£'Lv(0)l)'^^^ where the electron wave function 
is mostly asymmetri c, while the transverse velocity of the ato m 
reaches a maximum dVincke et al.lll992l:|Porekhinlll994[l998h . 

3.1.2. Continuum 

Wave functions of the continuum are computed using the same 
expansion, Eq. ( fT9] l. and system of equations in Eq. ( l20b . as for 
the bound states, but for a given energy E for every z-parity. The 
soluti on is based on a translation of the usual /^-matrix formalism 
(e.g.. lSeatonlll983b to the case of a strong magnetic field. Now 
\k) = \ns, E, K, +), where "+" reflects the symmetry condition 
g(z) - +gi-z)- Numbers n and s mark a selected open channel, 
defined for E^ s E - E^^ > by asymptotic conditions at z ^ 

sZ^I;k(z) ~ Sn„nSs„s sin[0„,(z)] + R„,,s„:„s cos[(f>„^^,^{z)], (22) 

where the pairs (ng, Sg), as well as (n, s), relate to the open chan- 
nels {E>E^,), 



siz) - KsZ + (m^e^/h^k„,) ln(^„sz) 



(23) 



is the z-dependent part of the phase of the wave function at 
z — > +00, and k„s - -y/2m,(£' - E^^)/h is the wave num- 
ber. For the closed channels, defined by the opposite inequality 
£';J- J > E, one should select g„^.s,.-K{z) — » at z — > oo. If /„ is 
the total number of open channels at given E, then the set of so- 
lutions, defined by Eqs. ( l20t and (1221) . constitute a complete set 
of /(, independent real basis functions. The quantities Rn',„s'„;n„s„ 
constitute the reactance matrix "R, which has dimensions /„ x /„. 
If the wave functions are normalized according to the condition 



JW J-L-/2 



dz\Mr^,zr^h 



then the reactance matrix satisfies the relation 

k , ,R , , — k R , , 



(24) 



(25) 



which differs from the usual symmetry relation (ISeatonlll983h . 

The representation with t] - \ must be used for continuum 
states, to ensure that the right-hand side (rh.s.) of Eq. (l20l l van- 
ishes at |z| — » oo, which is required by the asymptotic condition 
ofEq. (123. 

For a final state of a transition, one should use wave functions 
describing outgoing waves. The basis of outgoing waves with 
definite z-parity is defined by the asymptotic conditions 

8Z.,nsiz) ~ 6„„„6,„, e'^-(^) - 5,;,^^„, e-'^""-«, (26) 

where Sn„s„;ns are the elements of the scattering matrix S - 
i\ + VR){\ - VRy^ . The matrix S is unitary, but (again unlike 
conventional theory) asymmetrical. The basis of outgoing waves 
is obtained from the real basis by transformation 



(27) 



Here, pairs («, s) and («', s'), being respective parts of the com- 
posite quantum numbers k and /c', run over open channels, but 



in" , s") run over all (open and closed) channels. From the uni- 
tarity of the scattering matrix it follows that the wave func- 
tions that satisfy the asymptotic condition given by Eq. ( |26] | 
should be multiplied by a common factor {2L,)-^I^, to ensure 
the normalization expressed by Eq. (l24l) . As the initial state 
of a transition, one should use the basis of incoming waves. 

After the ortho-normalized outgoing waves have been con- 
structed for each z-parity, with symmetric and antisymmetric 
longitudinal coefficients ?°"'„ i-r^(z) - ±;?°"'„ ,^,^A-z) in 

'^n" s"\nsKE± '^n" s"\nsKE± ' 

expansion (T9[ . solutions for electron waves propagating at z — > 
±oo in a definite open channel (n, s) with a definite momentum 
hk - (signz) fik„s are given by the expansion in Eq. ( fT9] ) with 
coeflicients 



?n" .':";nskiz) — (.8 n" s" ■,nsE+^-^^ - 



' s"\nst- 



(28) 



where the sign -H or - represents electron escape in the positive 
or negative z direction, respectively, and we have suppressed K 
in the subscripts. Waves incoming from z — » ±00 with a definite 
momentum are given by the complex conjugate of Eq. ( |28] ). 

3.2. Adiabatic approximation 

In early works on the H atom in strong magnetic fields, 
a so-call ed adiabatic approximation was wid ely used (e.g., 
lGor'kov & Dzvaloshinskii 1968: Canuto & Ven tura 1977 and 
references therein), which neglects all terms but one in Eq. (fT9] l. 
i.e.. 



Mr) = '^ns{r±)g,{z). 



(29) 



This approximation reduces the system ( l20b to the single equa- 
tion with (n", s") - (n, s) and zero r.h.s. 

The accuracy of the adiabatic approximation for bound states 
can be assessed by comparing \eI\ with the distance between 
the neighboring Landau levels that are coupled by the rh.s. of 
Eq. (I2OI 1. For an atom at rest (K - 0), all the channel-coupling 
terms become zero for s" + s. In this case, the relevant Landau 
level distance is huj^, while the longitudinal energies of the states 
wifli V = ("tightly bound states") are I^U.f/O)! ~ 0.1-0.3 keV at 
B ~ 10"-10'"^ G , so that the adiabatic approximation is accurate 
to within a few percent or better It becomes still better for the 
"hydrogenlike states" with v > 0, which have \eI\ < 0.01 keV. 

From comparison of [^^^^.(O)! with hQ.c, one can conclude 
that the adiabatic approximation is generally inapplicable to a 
moving atom. However, the accuracy remains good for suffi- 
ciently slow atoms, that is when either <s: K^i or (provided 
that 77 = 1) » K„ (lPotekhinill994lll998h . Otherwise, since 
off'-diagonal effective potentials V„"s".n'.^^'(fB,z) in Eq. (|20] | de- 
crease at |z| — > CO more rapidly than diagonal ones, this approx- 
imation accurately reproduces wave functions tails at large |z|, 
provided that t] - I. 

For continuum states, the reactance and scattering matrices 
are diagonal in the adiabatic approximation, with a separate scat- 
tering coefficient S„s - S„s-„s for every open channel. 

3.3. Born approximation 

In the Born approximation, the potential V in a Hamiltonian 
H - Hi) + V, which acts on particles in the continuum states, 
is treated as a small perturbation. We define i^**^* to be the non- 
perturbed function, which satisfies the equation //qi/'**'' = £1/'*"*. 
Then from the Schrodinger equation Hiff - Eif/, one obtains the 
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continuum wave function in the first Born approximation in the 
form i/r = i/r**^* + t/r*'*, where i/'*'* is determined by the equation 
(Ho - E)tJ/^^> = -yi/r(">. 

Since we consider the continuum states corresponding to 
definite Landau numbers (n, A^) at |z| — » oo (Sect. l3.lT2] i. the zero- 
order wave function i/'*"* is given by the adiabatic approximation 
with giiiz) replaced by plane waves. 

3.3.1 . Two forms of solution 

We now consider continuum states. We first choose the non- 
perturbed wave function in the representation where the z- 
projections of angular momentum operators, have definite 
values -hse for the electron and hsp for the proton. Then 

'^!,°.L-t,„w,.pA^('-e' '•p) = 'i'n,...kJr.)'i'N,s,.kJrp), (30) 
and i^*'* is governed by the equation 

(i/<" + i/« - e) .A<"(re, rp) = (eVrep) <A*°'(re, r,). (31) 

Using an expansion of t/r*'* over the complete set of 
k' N' s' k[ ' obtain in the standard way 

^ n' ,s' ,N' ,s" ,k'. 



where k, = {k,f.M - k~pm)/mu - k,s - -k~p (since P, = 0) and 
dze'^<^ drldr, 

KA^l) ^N'."(r±) ^nsMl) %s(''-l) 



(|rl-rj2+z2) 



1/2 



In the limit L; — > oo, we replace by (L^jln) dk,. 

fPotekhin & Chabrier ( 2003) obtained a simpler solution, 
based on the representation of quantum states with definite K. 
In this case, there are no separate quantum numbers and s^. 
After applying the transformation in Eq. ( fT3]) . ^'^"^ is given by 
Eq. ( |29] l with gxiz) - e^^'*^--/ ^|Ll. Using Fourier transform 



1 rL~J2 
V^z J-LJ2 



(34) 



we obtain from Eq. ( l20l l in the first Born approximation 

<k.(r^, k) = ijjfKr^, k) + 4^['\r^, k) (35) 
with 

~((\\ sin((^ - ^.)L,/2) ^ ^ ^ 



{k-kdLJ2 
- (271/ L,) 6(k - k,) ^ns(r±) at ^ oo, 



V:'J (r,,k) 



-L:'V':/(r,,k-L) 
n\e-e)l(2m,) + E^,,-E^, 

e"'*' V„„,'.'(rc,z)dz, 

OO 



(36) 
(37) 

(38) 
(39) 



where V"/ can be presented as a single int egral of a combination 
of ele mentary functions (Appendix B of iPotekhin & Chabrieii 
I2003h . 



3.3.2. Approximation of infinite proton mass 

The neglect of the proton motion is equivalent to the assumption 
that M — > oo. In this approximation, i/'*"* depends only on in 
Eq. (l30b without ^N,s^,Lf.i^v^ on the r.h.s. Then Eq. ( l32b simpli- 
fies to 



<A^"('-e) 



n^X. E^,+(hk';)^/2m-Ej_,^ ' 



(40) 



Taking into account the definition in Eq. (|2TJ, we see that this 
solution is identical to the solution provided by Eqs. (l37])-([39l) 
in the particular case where = 0, after the obvious replacement 
of E^^ by E^ and m, by m. The zero value of rc naturally reflects 
the condition Vdiift = 0. 

4. Electron-proton photoabsorption 

4.1. General expressions 

The general nonrelativistic formula for the differential cross- 
section of absorption of radiation by a quantum-mechanical sys- 
tem is (e.g., Armstrong & Nicholls 1972) 



(32) do- : 



4n^ 



\if\e-j.fi\i}f 6(Ef-Ei-hio)dvf, 



(41) 



where |/) and |/) are the initial and final states of the system, dvf 
is the density of final states, ha> is the photon energy, e is the 
polarization vector, jeff - e'* ""j, j is the electric current operator, 
and q is the photon wave number In our case. 



(33) yeff = e(e"''-'re-e'*'-prp). 



(42) 



where the velocity operators and are given by Eq. (O. 

Equation ( |42] ) does not yet include either the photon inter- 
action with electron and proton magnetic moments ji^ or ft^. 
For transitions without spin-flip, the latter interaction can be 
taken into account by adding to the e ■ / "etr o perator the term 
jspin = -i(9 xe) ■ (fi^ + pip) (cf. lKopidakis^traLHl996h . whereas 
ope rators (g xe) X j j.„ ^ are responsible for spin-flip transitions 
rcf. lWunneretal] IT983). 

We consider the representation where and Sp are definite 
in the initial and final states. For an initial state with fixed n,-, s^j, 
Ni, Spj, and L = k, in Eq. (l32T l. and for a final state with either 
a fixed z-parity or a fixed sign of k- = kf, we have in Eq. (HTt 
dvf = (L^/2n)(mt/h?\kf\)dEf. Therefore, the cross-section of 
photoabsorption for a pure initial quantum state \i) is 



ri(co) =2 



" \{f\J\i)\\ 



h?-\kf\a>c 

where J = e- jeft + Jspin, 

ffk^ 



(43) 



E^ 



2m, 



= E-^ 



2ot, 



-I- ha), 



(44) 



the sum is performed over those n j and Nf which are permitted 
by Eq. ( l44l i. and "+" means the sum over the z-parity of the final 
state (in the case where the parity is definite) or the sum over the 
signs of kf (in the case where k^ in the final state is definite). 

In the alternative representation with definite cartesian com- 
ponents of pseudomomentum K, using the transformation in 
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Eq. ( fTlT l. one can express the cross-section in terms of the inter- 
action matrix element between the initial and final inte rnal states 
of the electron-proton system ( Bezchastnov & Potekhinlll994l) . 
The result has the same form as Eq. (HTt . but now dvf is the den- 
sity of final states at fixed Kj - Kj + hq, initial and final states 
are described by wave functions tfru, and the effective current op- 
erator in the conventional representation with tj - (r - r^p) is where 
given by 



By substituting Eqs. ( |46] | and ( |47] ). the sum in Eq. dSOb reduces 
to 



Qr--1 n's' 



JeS 



M \ln K hq 
e expji — q r\\- + - — + — 



•f n' s' 



ihel yjn' + 1 



^ n'N' 



.N'-\ 



,exp,-,-,.r|^H_ 



K _ hq_ 

hmii 2M 



(45) 



n s 



■L- 



n'N' 



M 



-n',N' + \ 
n'N' 



where operator n is defined by Eq. ( fTTb and H = p- (e/lc) Bxr. 
The transverse cychc components of operators n and II act on 
the Landau states as 



ih 



(53) 

(54) 
(55) 
(56) 
(57) 



(46) 
(47) 



where = n + j ± j and N"^ = N + j + j - n"^ + s. 

Changes in and induce transformati ons of operator 
Jeff, studied bv iBezchastnov & PotekhinI (|1994|) . In the particu- 



and = n + i is the proton Landau number In the first Born 
approximation (Sect. 13.3] ). 



{f\J\i)^ <«A/ * I e ■ Jeff I ifrf^) + {ilff I e ■ Jeff I .A*")- 



(58) 



In the representation where ^e and Sp are definite, using Eq. ( l30b 

for 1^**^^ and Eq. (|32] | for i^*'', and taking into account the relations 

1 u r 1- .1- ■ -^^ 1 J £ 1 4. 4. 4 4- in Eq. dgll, we can derive the explicit expression for the matrix 

lar case where for both mitial and final states, the representation gjgjj^jjfjjj g; d^iK of f f 

with T] - 1 {vb - rc,r - rep + r^) is used, their result reads ^' ^ — ' 



— - exp I iq ■ 



' ep 

r^+ra M . 

— ^ — + — ^qzZ 
1 mil 



+ exp -iq 



m 



n hq 

m 2m 



n 



+1 



hq 



i?-^^ 77-:T77 ■ (48) 



M 2M 



In this representation, instead of Eq. ( l43T l. we have 



where 

;(0) _ 



0-, 



Z2nL.mt ^ 



(49) J 



:(±1) 



L~ m, CO 
ie 



o=-l 



(59) 



(60) 



1 



Ljam [m(Lo± Wc) 



„± TJi7"f^l'^e.f'fl,Nf.Spj 



(Ak) 



where the sum is performed over those «/ and Sf that are per- 
mitted by Eq. (l44l i. For the solution described in Sect. 13. II the 
matrix element in Eq. ( |49] l becomes 



1 



{f\J\i} 



Z [%Ts''-,.M 

x{n"s"\J\n's'}^g'^^,.Jz)dz. 



M(aj + ilc) 

n/,s„,A',±l,.SpjT 



(61) 



and Ak = kf - ki. 

In the representation where cartesian components of K have 
definite values, using Eqs. (l37])-([39ll, ( |46] |. and ( |47]) . one can de- 



(50) 

Using Eqs. (l46l l and (l47l l. we can express the transverse ma- 
fiix elements {n" s" \J'\n's')i_ in terms of Laguerre functions. 

Hence, Eq. ^ presents a s um of overlap integrals ov er z. For rive the mafi-ix element in Eq. (gUl in the form 
instance, Eqs . ( A7 )-( A 1 2) of lPotekhin & Pavlov! (Il997h provide + 1 

an explicit expression in terms of this overlap integrals for the 
matrix elements {f\M \ i) of the operator M - (hj e^) J' \n the 
approximation where small terms ~ OiimjM) q) are neglected, 
but separate terms ~ OimjM) and ~ 0(q) are retained. 



4.2. Dipole and Born approximations 

Hereafter, we use the dipole approximation (q — > 0). Then jspin 
vanishes, and the total effective current in Eq. (l42l) reduces to 

Jeff = e(re - j-p), (51) 

while the transformed effective current in Eq. ( l48b becomes 

Jeff = e(;r/m + n/M). (52) 



where 

J(0) ^ 
7(±1) - 



L~JJL(x) 



AkVZY(r„Ak), 



(62) 



(63) 



le 



[V^V-C'-^'"(..,A^)-V^V:::i;,^,,(r.A^) 



m (W + We) 



■N-VZ:'^\r,,Ak)- 



wv„:::,,(rc,Ak) 



M{aj + Qc) 



(64) 
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Substituting Eqs. (|62]|-(|64]| into Eq. ( l49b , assuming Maxwell 
distribution of^fcj, and taking the average over the initial states, 
we obtain JPotekhin & Chabrie J2063t IPotekhTn & Laili2007b 



+1 



(65) 



a=-l 



where and are the electron and proton number fractions at 
the Landau levels n and A^, 



cr, 



(a) 



,2 „(a) 



O) V, 



n,N; n',N' 



(66) 



is the partial cross-section for transitions between the specified 
electron and proton Landau levels for polarization a. 



^,M;„',W'(^) 




A 



{a) 

n,N;n',N' 



mT fui) 

is the effective partial collision frequency, 



(67) 



n,N;n',N' 



(/?,, oj/co, 



2 Jo W 



X ((«' + M)2i«l w%.„, „,iu^) + iu' - «)2i«l w%,, ,,.(u^)) (68) 
is a partial Coulomb logarithm, and 

df 

(f + Mi/2)2 



1 r" 

>^11;„',A^'(«±) = 2 Jo 



it). 



where y6» = hojt/T - fieB/(m,cT), u± - \u + u'\, 9{u' ) is the 
Heaviside step function, and 



,7 o 2m.t , 2m, , 2a» 

M m ojt 



(70) 



Since A*^'^ = A*"'', two different Coulomb logarithms A*"' and 
A*'''' describe all three basic polarizations. 

Terms that are proportional to e^e*, with a ff' are absent 
in Eq. (l66T l. because, for every pair of pure quantum states |/) and 
I/), only one of the three basic polarizations provides a non-zero 
tra nsition matrix element in the dipole approximation. 

iPotekhin & Lail ( l2007b mentioned that Debye screening 
might be taken into account by using M± - [(M+M')^+(flm^D)^]'^^ 
as the arguments of w'"]^. ^, in Eq. (l68T l. ko being the inverse 
screening length. However, ISawved (l2Q07h . following iBekefil 
showed that scattering off a Debye potential is not a 
valid description of the screening correction for photoabsorp- 
tion; instead, the integrand in Eq. ( |69] | should be multiplied by 
(t + + fl^fc2^/2)/(f + + a^k^), where k^^ is the electron 
contribution to the squared Debye wave number k^. 

4.3. Damping factor 



Equation (166b gives divergent results at w ^ for a - -1 and 
at w ^ Qc for a - +1, because it ignores damping effects due 
to the finite lifetimes of the initial and final states of the transi- 
tion. A conventional way of including these effects consists of 
adding a damping factor to the d enominator in E q. (|66l), whic h 
results in Lorentz profiles (e.g.. [Armstrong &~Nichollsill972h . 
The damping factor can be traced back to the accurate treat- 
ment of the complex dielectric tensor of the classical magnetized 
plasma CGinzburail970.) . This treatment allows one to express 
the complex dielectric tensor in terms of the effective collision 



frequencies related to different types of collisions in the plasma. 
Imaginary parts of the refraction indexes, calculated from the 
complex dielectric tensor, provide complicated expressions for 
the free-free photoabsorption cross sections cr^ for the basic po- 
larizations a = 0, +L Based on the assumption that the effective 
collision frequencies are small compared to oj, the latter expres- 
sions greatly simplify and reduce to ( P otekhin & Chabrier 20031) 



where 



a? vff(w) 



mc (a> + atOc)^ (a> — aD.c)^ + lo^ v^(tj) ' 



(71) 



(72) 



Ve and Vp being the efl^ective damping factors for protons and 
electrons, respectively, not related to the electron-proton col- 
lisions. In gener al, Ve and Vp may also depend on a and a>. 
jjinzbur^ d 19701) considers Ve and Vp fo r comsions of elec- 
trons a nd protons with molecules, whereas fPotekhin & Chabrial 
(l2003h take into account damping factors due to both the scat- 
tering of light by free electrons and protons and proton-proton 
collisions. The derivation of Eq. ( l72l i from the complex dielec- 
tric tensor of the plasma assumes that Vg <k oJc, Vo ^ Wc, and 
Vp <K Q.C. 

Although the general expressions given in Eqs. ( fTTT l. (|72] ) can 
be established in frames of the classical theory, accurate values 
of the effective frequencies are provided by quantum mechanics. 



(59) In our case. 




(73) 



where vj^"^.^, ^,(w) is provided by Eqs. dFTb-dTOl). In the second 
equaUty, A^ and are, by definition. Coulomb logarithms for 
a = and a = ±L Parallel and transverse Gaunt factors (e.g., 
lMeszaroslll99Z) equal ( VI/tt) A^^ and ( VS/tt) A^, respectively. 

Since diff'erent quantum transitions contribute to the cy- 
clotron resonance at the same frequency (Qc or Wc, depending on 
a), their quantum amplitudes are coherent. Therefore it is impor- 
tant that the same damping factor Vq, be used in all the transitions 
(cf. the discussion of radiative cascades in quantum oscillator by 
Cohen-Tannoudii et al. 1998). Moreover, the same v„ given by 
Eq. ( |72] | should be used for the absorption and scattering pro- 
cesses. This ensures that the cyclotron cross-section, being in- 
tegrated across the resonance, pro vides the corre ct value of the 
cyclotron oscillator strength ('e.g.. lVenturalll979 i). otherwise the 
equivalent width of the cyclotron line would be overestimated. 

In the electron resonance region, where |w - wd co^ and 
a - -I, one can neglect Q.c/a>, because it is much smaller than 1, 
and the term that contains Vp, because it is small compared to the 
other terms. The result coincides with the conventional expres- 
sion for the electron free-free cross-section without allowance 
for proton motion with v_i - v^^ -H Vg. In the proton resonance 
region, where |w - Qcl ^ ^^c and a = +1, the denominator in 



Eq. dTTTi becomes (a>+ojr)^(tD-Q.r)^ +tdP-vl ^ ul [(w-Qc)^H-Vp], 
where Vp - (m/M) v|^(Qc)- In this approximation, Eq. (fTTT l be- 
comes formally equivalent t o a simple one-par ticle cyclotron 
cross-section (cf. Eq. 14 of iPavlov et al.l 1 19951 or Eq. 47 of 
ISawveiil2007l) . apart from a difference in notations and the dif- 
ference in a"" (the latter being discussed in Sect. l5.2b . 

The treatment that leads to Eq. ( |72] | predicts a small shift in 
the position of the resonance due to the damping. This shift is 
unimportant for applications and therefore neglected in Eq. ( fTTT l. 
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5. Cyclotron harmonics 

In addition to the fundamental cyclotron resonances, the quan- 
tum treatment of the free-free absorption identifies electron and 
proton cyclotron harmonics at integer multiples of Uc and Q.^, 
respectively. They appear because of the increase in the partial 
Coulomb logarithms Ajf^.^, ^,(y6,, w/w,) at w ^ Oc(A^' - A^) + 
(jJcin' — n). Thus, /th electron cyclotron harmonics (in addition 
to the fundamental at w = Wc) arises at w = {I + 1)0)^ due to 
the terms with n' - n - I + \, and each Zth proton cyclotron har- 
monics (additional to the fundamental at w = Qc) is formed by 
the terms with A^' - = / + 1 in Eq. ( l73T l. Unlike the classical 
electron and proton cyclotron resonances, the quantum peaks of 
A*"^ contribute to cTaico) at any polarization and are the same for 
a - +\ and - 1 . 

The relative strengths of the harmonics depend on the dis- 
tribution numbers /„'' and In this paper, we assume local 
thermodynamic equilibrium (LTE) and thus use the Boltzmann 
distributions, as in most of the previous papers (but see 
[Nagel & Ventura 1983 and Potekhin & Lai 2007 for non-LTE 
effects on the electron and proton cyclotron radiation rates, re- 
spectively). 

We calculate free-free cross-sections in magnetized neutron- 
star atmospheres using Eqs. (ISSll-dTSll. Examples of opaci- 
ties and/or spectra calculated with the use of these cross- 
sections can be four id, e.g., in [Po tekhin & Chabrier (2003, 
12004. iPotekhin et all (l2004 . lHo etaLi (,2008.) . ,Suleimanov et ah 
(120091 12010ar In previous studies, various additional simpli- 
fications have been made in addition to the nonrelativistic, 
dipole, first Born approximations described above for the free- 
free cross-sections. Below we assess the applicability ranges of 
these simplifications by comparing with our more accurate re- 
sults. 



5. 1. Electron and muon cyclotron harmonics 
5.1.1. Fixed scattering potential 

In early works (e.g., Meszaroslll992l and references therein), 
free-free (or bremsstrahlung) processes were treated assuming 
scattering off" a fixed Coulomb center, which is equivalent to the 
approximation of M — » oo, described in Sect. l3.3.2l In this ap- 
proximation, one can set Oc = and explicitly perform the sum- 
mation over A^' in Eq. (l65l l using the identity Yjn'=o ^n'n^^^ ~ ^■ 
Taking damping (Sect. l4.3l l into account, we obtain 



Ane^ 



mc (oj + aojc) + (ve + vjj)- 



S\2 



ff 4 / 27r He k 
3 V ml fujj 

(^e, Ulu.) - 5 Z Z r ^ 

n n' 



x((«' + «)2Nwl'^>,(«.)- 



f I"! df 



(f + «i/2)2 



(74) 
(75) 

(76) 
(77) 



where [5^ - ha>c/T and u - + 2 (n — n') + IojIo)^. Assuming 
Boltzmann distribution (f^/f^ - 2e^"^' at « > 1, where the fac- 
tor 2 takes account of the elect ron spin degeneracy), one can 
reduce th is result to Eq. (27) of | Pavlov & Panovl d 19761) (as cor- 
rected bv lPotekhin & Chabrieill2003h . 
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Fig. 1. Transverse Coulomb logarithm as a function of a>/cjc 
at fi,(jjc/r = 5 for different approximations: the model of fixed 
Coulomb potential (dotted line), approximate account of proton 
recoil (solid line), adiabatic approximation (long-dashed line), 
and the first post-adiabatic approximation (short dash - long 
dash). The divergent peaks are trimmed at \a) - (l + = 
lO^^'wc (I - 0,1,2,...). For comparison, the nonmagnetic 
Coulomb logarithm (short dashes; in this case the horizontal 
axis displays ha)/5T) and the model with approximate account 
of proton recoil in the muonic atom ju p (dot-dashed line) are 
shown. 
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Fig. 2. The same as in Fig.[T|but for the longitudinal Coulomb 
logarithm. In this case, lines with the approximate account of 
proton recoil almost coincide with the line corresponding to the 
fixed-potential approximation. 
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This approximation was used in all models of the spectra of 
strongly magnetized neutron^stars until the beginning of the 21st 
century (e.g.. , Pavlov etal.1 19951 : IZane et al. 2000, and references 
therein). It is validated by the large value of the mass ratio M/m. 
In addition, it requires that oj » fic, as seen directly from the 
comparison of Eq. ( l74l i with Eq. dTTT l. 

5.1.2. Approximate account of proton recoil 

iPavlov & Panovl d 19761) proposed an approximate treatment of 
proton recoil, which assumes that j6e «: M/m and does not 
take into account Landau quantization of proton motion. In 
Fig. [U the dotted line shows the perpendicular Coulomb loga- 
rithm A*[ calculated according to Eqs. (I74ll-(f77]l. while the solid 
line takes the approximate account of proton recoil. As an ex- 
ample, we show the case where [5^-5. The familiar nonmag- 
netic Co ulomb logarithm in the first Born approximation (e.g., 
iBethe & Salpeter 1957) is shown by the short-dashed line, as- 
suming the same ftw/T along the horizontal axis as for the other 
curves (ha>/T = /Ssco/coc)- 

To enhance the difference caused by the recoil, we replace 
the electron by the muon // . All the above formulae and discus- 
sion remain unchanged, but now the mass ratio is M/m - 8.88. 
The result of the approximate treatment of the recoil is shown by 
the dot-dashed line. 

In Fig.|2l the same approximations are shown for A^^. In this 
case, the lines related to the cyclotron harmonics are smoothed, 
because the factors (m' ± u)^ quench the near-threshold growth 
of the integrand in Eq. ( |76] |. The same smoothing results in the 
infinite proton mass approximation being even more applicable 
(under the necessary condition co » Q.^)' the dotted, solid, and 
dot-dashed lines almost coincide in Fig.|2] 

5.1 .3. Adiabatic and post-adiabatic approximations 

Several authors (IVirtamo & J auholl975l:lNag el & Venturall983t 
eszarosl[l992l) used the adiabatic approximation not only for 
the unperturbed wave function but also for i//^^\ This 
was done in addition to assuming the infinite proton mass 
(Sect. l5.lTn i. In other words, they kept only one (n, s) term in 
the sum given by Eq. ( l40l l. The result is shown in Figs. [T] and |2] 
by long-dashed lines. We see that this approximation works well 
at w_«_Wc, but becomes inaccurate at > Wc- 

ISaw ver (2007) analyzed the photoabsorption problem by us- 
ing the method of field theory. In the region Q.^ <^ cj < I .Sate, he 
considered account two electron Landau levels n = and 1 and 
applied a perturbation theory assuming the parameter s = e"^'^^ 
to be small. His result is identical to the results discussed in 
Sect. 15 . 1 . TI expanded in powers of s, which we can write as 

+ ——[Ko(xi) + Ko(x^i))j. (78) 

Here and in the next equation, Ky(x„) are modified Bessel func- 
tions, and x„ = \ h<x)/T + nPe\ -s/0 .25 + y//3e. Equation dTST l differs 
from Eq. (28) of ISawve^ (l2007h in two respects: first, we have 
restored the factor 2 at erP'l^^ and second, we have dropped a 
term proportional to e"^% because it is of the same order as 
the contribution from the level n - 2, and therefore should be 
treated together with the latter contribution in the next order of 
the perturbation theory. 




0.5 1 1.5 2 2.5 



Fig. 3. Transverse Coulomb logarithm as function of w/wc at 
hcOc/T = 5 in the model of fixed Coulomb potential (dotted 
line) and the approximate account of proton recoil (dot-dashed 
line), compared to the nonmagnetic Coulomb logarithm (dashed 
line) and the accurate calculation for the muonic atom p (solid 
line). 

In the same way, we obtain 

A!f = -e'"-/^ r ^ lx,Ki(xo) 

" 4 Jo (l+y)if^J4 + y)\ 

+ j-^{xiKi(xi) + x^iKi(x^i))j- (79) 

Equations ( l78l l and (|79]l can be obtained by the first iteration in 
the perturbation expansion, starting from the adiabatic aproxi- 
mation. 

5.2. Proton cyclotron harmonics 

Proton cyclotron harmonics in the photoabsorption coefficients 
at w = (/ -H 1) Qc are superimposed on the peaks related to the 
electron cyclotron harmonics. However, for the H atom the two 
series of harmonics are separated because of the large value of 
M/m = 1836.1. To observe the superimposition and the quali- 
tative differences of various approximations, it is instructive to 
consider, in place of the H atom, the muonic atom (the fi^p sys- 
tem), which has a smaller mass ratio M/m - 8.88. The trans- 
verse Coulomb logarithm of photoabsorption by such sys- 
tem is shown in Fig. [3] The solid line displays the result of a 
calculation made according to Sect.|4] The other lines, as well as 
in Fig. [T] show the results of different approximations: a fixed 
Coulomb center (Sect. l5.lTTl dotted line), the approximate ac- 
count of proton recoil (Sect. 15.1^ dot-dashed line), and a non- 
magnetic Coulomb logarithm (dashes). 

The smaller peaks in the solid curve correspond to the proton 
cyclotron harmonics. They are superimposed on the large-scale 
oscillations, which correspond to the muon cyclotron harmonics. 
Although the approximate recoil treatment (dot-dashed line) im- 
proves the agreement with the exact calculation compared to the 
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Fig. 4. Transverse Coulomb logarithm as function of a>/Q.c at 
hQc/T = 5. The accurate calculation (solid line) for the systems 
jU^p (lower lines) and ep (upper lines) is compared to the approx- 
imation of a fixed Coulomb potential for the electron or muon 
scattering (dashed lines) or for the proton scattering (dotted line), 
and to the first post-adiabatic approximation (short dash - long 
dash). 



infinite proton mass model (dotted line), both that approximate 
models that neglect proton Landau quantization differ signifi- 
cantly from the precise result. 

In Fig. H] we compare the proton cyclotron harmonics for 
different relative masses of the positive and negative particles. 
Here the proton cyclotron parameter is fixed to /5p - Ml^T - 5, 
and the horizontal axis displays the ratio co/Q.^:. The solid lines 
show the transverse Coulomb logarithm for the muonic atom 
(the lower curve) and the H atom (the upper curve). The dashed 
lines show A|[ calculated for the same TuDdT - (3pM/m and 
the same lj/cOc - (/«/M)(y/Oc as the solid curves, but for the 
approximation of a fixed Coulomb potential in the electron or 
muon scattering. By comparison, the dotted line shows cal- 
culated for proton scattering off a fixed Coulomb center, which 
can be regarded as a model where M/m — > 0. We see that the ap- 
proximate models are unable to reproduce the proton cyclotron 
features correctly. It is also noteworthy that the larger the ratio 
M/m, the smaller the proton cyclotron peaks. In addition, the 
cyclotron resonance strength decreases with increasing harmon- 
ics number /. These properties of the cyclotron harmonics allow 
us to conclude that the solid lines in Figs. [T]and|2]are precise 
(proton cyclotron harmonics are negligible on their scale). 

In the ea rly models of mag netized neutron star atmo- 
spheres (e.g., Pavlov & Shibanov| |1978t IShibanov et al] 119921: 
IShibanov & ZavUn J^5), the authors considered moderate 
magnetic fields B ~ 10"-10'^G, where the proton Landau 
quantization is unimportant. More recently, observational evi- 
dence has accumulated that some of the isolated neutron stars 
are probably magnetar s, which have fi elds of B ~ 10'"* G (see, 
e.g., the review by iMereghettil |2008| and references therein). 
According to Eq. the proton cyclotron lines of magnetars 



are in an observationally accessible spectral range, which has en- 
couraged theoretical modeling of these features. In t he absence 
of ari accur a te quantuni treatment, several auth ors dZane et al.l 
I2000L 1200 U lOzeil 120011: IHo & Laill200lL l2003h e mployed the 
scaling previously suggested for this purpose by iPavlov et aH 
([T995r) . according to which the free-free cross-section for pro- 
tons equals (m/M)^ cr_Q,(wM/m), where cr„(w) is given by 
Eq. (f74l l. The latter equation differs remarkably from the cor- 
rect expression in Eq. (fTTT i. At photon frequencies cj < Oc, the 
difference roughly amounts to a factor of {co/D.^}^. 

In addition, the Coulomb logarithm that determines can- 
not be obtained from this scaling. An example is shown in Fig.|4] 
where the dotted line corresponding to the fixed-potential model 
is compared with the accurate calculations displayed by the solid 
lines. We see that the fixed-potential model strongly overesti- 
mates the strength of the proton cyclotron harmonics. The origin 
of the discrepancy is clear: while considering a collision of a 
proton with an electron, one cannot assume the electron to be a 
nonmoving particle. 

ISawver (2007) employed a representation with definite ie 
and Sp and analyzed the first proton-cyclotron peak of A^, in a 
way similar to his analysis of the first electron cyclotron peak 
(see Sect. l5.1.3T l, by taking into account the ground electron 
Landau level n - and two proton Landau levels, N - and 1 . 
The result (his Eq. 30) is quite accurate close to the fundamental 
cyclotron frequency, as we illustrate wuth the lines of alternating 
short and long dashes in Fig. |4] In the case of hydrogen (higher 
M/m), it almost coincides with the accurate result (solid line) 
at (L) < 1 .5 Qc and with the result obtained by neglecting the 
Landau quantization of protons (dashed fine) at higher w values. 

6. Discussion 

6.1. Corrections beyond Born approximation 

The formulae presented in Sects. [TT] and 14. 1 1 in principle allow 
one to perform an accurate calculation of photoabsorption rates 
in the electron-proton system in an arbitrary magnetic field, tak- 
ing into account the effects of Landau quantization of the elec- 
tron and proton motion across the field and the transverse mo- 
tion of the center of mass. For bound-free ab sorption, this cal- 
culation was presented bv iPotekhin & Pavlov! (fl997.) . For free- 
free processes, we apply the first Born approximation and the 
dipole approximation. We plan to perform calculations of the 
free-free opacities beyond Born approximation in future work. 
An appr oximate estimate of th e non-Born corrections can be ob- 
tained ( Potekhin & Laill2007h by introducing correction factors 
(1 + r;'M"2)-'/^(l -I- 7:\uT^)'^'^ into the integral of Eq. dSH), 
where y, = {m/m»fy and y = h^B/(m^ce^) = 425.44 B^. The 
accuracy of the approximation is ensured by the smallness of 
r;'^^ ~ 0.05 5^2^^ additional condition T » e^m jh^, 

which is the usual applicability condition for a Born approxima- 
tion without a magnetic field. 

We have checked that these corrections are sufficiently small 
for the electron cyclotron harmonics at B ~ 10" G (relevant 
to CCOs) and negligible for the proton cyclotron harmonics at 
B > 10'^ G (relevant to XDINSs). 

6.2. Importance of bound States 

Free-free absorption contributes only a part of the total opaci- 
ties in the atmospheres of neutron stars. A second constituent 
is the familiar sc attering, and a third the absorption by boun d 
species (see, e.g.. ICanuto & Venturalll977l : iPavlov et al.lll995h . 
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IogB=13.7, T=120 eV, p=10 g cm-3 




0.2 0.4 0.6 0.8 1 2 

hcj (keV) 

Fig. 5. Opacities for the two normal electromagnetic waves 
propagating at the angle 9b - 45° to the magnetic field direction 
in a hydrogen atmosphere of a neutron star with Z? = 5xlO'^G 
and T = 120 eV at density p = 10 g cm"^ (which is in the mid- 
dle of the photosphere at these B and T). The results are shown 
for fully ionized (dotted lines) and partially ionized (solid and 
dot-dashed lines) plasma models. In the latter model, the non- 
ionized atomic fraction equals 0.0066. The solid line shows the 
opacity obtained with the accurate calculation of the free-free 
Coulomb logarithm, and the dot-dashed line demonstrates the re- 
sult of the approximate treatment that corresponds to the dashed 
line in Fig.|4] 



It was realized long ago (iRudermanlllQ?! *) that in strong mag- 
netic fields the increase in the binding energies of atoms and 
molecules can lead to their non-negligible abundance even in 
hot atmospheres. With increasing B, the binding energies and 
abundances of bound species in creas e at any fi xed density p and 
temperature T (IPotekhin et al.lll999l;lLaill200lh . so that even the 
lightest of the atoms, hydrogen, provides a noticeable contribu- 
tion to the opacities at the temperatures of interest, if the mag- 
netic field is strong enough. Even a small neutral fraction can 
be important, because the bound-bound and bound-free cross- 
sections are large close to certain characteristic spectral energies. 

For electron cyclotron harmonics to appear at fi-w < 1 keV, 
we should ensure that B < 10" G according to Eq. ([1]). At these 
relatively weak magnetic fields and the characteristic tempera- 
ture r°° > 100 eV, the assumption of full ionization may be 
acceptable. However, at B > 10'^ G, which is required for ion 
cyclotron harmonics, the situation is different. An illustration is 
given in Fig.|5] The solid curves show true absorption opacities 
for two normal electromagnetic waves propagating at the angle 
45° to the magnetic field lines at B = 5 x 10'^ G and T = 120 eV. 
The upper and lower curves correspond to the ordinary and ex- 
traordinary waves, respectively. The density in this example is 
chosen to be p = 10 g cm"-*, which is a typical atmosphere den- 
sity at B = 5 X 10'^^ G and T^^ = 100 eV (at this density the 
thermodynamic temperature T approximately equals the effec- 
tive temperature T^s). According to our ionization equilibrium 



model frotekhin et"anil999l) . at these B, T, and p values, 0.66% 
of protons in the plasma are comprised in the ground-state H 
atoms that are not too strongly perturbed by plasma microfields 
so that they contribute to the bound-bound and bound-free opac- 
ities (the "optical" atomic fraction), and only 0.1% of protons 
are in excited bound states. Even though the ground-state atomic 
fraction is small, it is not negligible. In Fig.|5] at hco > 0.4 keV, 
the opacities in two normal modes, calculated with allowance 
for partial ionization (solid and dot-dashed curves), are signif- 
icantly higher and have more characteristic features than the 
opacity calculated under the assumption of complete ionization 
(dotted lines). In particular, the broad feature on the lower curve 
near 0.4 keV is produced by the principal bound-bound transi- 
tion between the two lowest bound states (sj = —>.?/ = 1), 
and the increased value of the opacity at higher energies hcj is 
due to the transitions to other bound and free quantum states. 
The wavy shape of the lower solid curve (for the extraordinary 
mode) at hco > 0.7 keV is explained by bound-free transitions 
to different open channels, each having its own threshold en- 
ergy. All the bound-bound absorption features and photoioniza- 
tion thresholds are strongly broadened by the effects of atomic 
mo tion across the magnetic field lines ("magnetic broadening", 
see lPotekhin & Pavlovll 19971 and references therein). 

In the insets, we zoom in on the regions of the first and sec- 
ond proton cyclotron harmonics. Both of them are visible, but 
negligible compared to the effect of partial ionization on the 
opacities. 



6.3. Other possibilities for CCOs and XDINSs 

Apart from the cyclotron harmonics, a number of alternative 
explanations of the observed absorption features in CCOs and 
XDINSs have been suggested in the literature. 

Im ori & Hoi (12007') constructed models of strongly magne- 
tized neutron star atmospheres with mid-Z elements and com- 
pared them to the observed spectra of the neutron stars IE 
1207.4-5209 and RX J1605.3-H3249. They demonstrated that the 
positions and relative strengths of the strongest absorption fea- 
tures in these neutron stars are in good agreement with a model 
of a strongly ionized oxygen atmosphere with B - 10'^ G and 
B - 10'^ G, respectively. This explanation seems promising, but 
unsolved problems remain: the effects of motion across the field 
have been treated approximately, based on the assumption that 
they are small, and detailed fits to the observed spectra have not 
yet been presented. 

Among other hypotheses about the nature of the ab- 
sorption features, there was a suggestion that they could 
be due to bound-bound transitions in exotic molecular ions 
dTurbiner & Lopez Vievrall2006i) . However, our estimates show 
that the abundance of these ions in a neutron star atmosphere 
would be negligible compared with the abundance of H atoms. 
Suleimanov et al. (2009) proposed a "sandwich" model atmo- 
sphere of finite depth, composed of a helium slab above a con- 
densed surface and beneath hydrogen, and demonstrated that this 
model can produce two or three absorption features in the range 
of h(jL> ~ 0.2-1 keV at B ~ 10'"* G, although a detailed compar- 
ison with observed spectra was not performed. One cannot also 
rule out that some absorption lines originate in a cloud near a 
neutron star, rather than in the atmosphere (see lHambarvan et al.l 
2009). 
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7. Summary 

We have considered the basic methods for calculation of free- 
free opacities of a magnetized hydrogen plasma. Our emphasis 
has been on the case where not only electron, but also proton 
motion across the magnetic field is quantized by the Landau 
states. We have derived general formulae for the photoabsorp- 
tion rates and considered in detail the dipole, first Born approxi- 
mation. We have presented numerical examples, compared them 
with the results of previously used simplified models, and an- 
alyzed the physical assumptions behind the different simplifica- 
tions and conditions of their applicability. We have demonstrated 
that the proton cyclotron harmonics at a given value of the pa- 
rameter ySp = fiQc/r are much weaker than the respective elec- 
tron cyclotron harmonics at the same value of jSg - hcjc/T, and 
explained this difference in terms the large (nonperturbative) ef- 
fects of proton motion in the case of proton cyclotron harmonics, 
in contrast to the case of the electron cyclotron harmonics. 
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